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Abstract 

The goal of this paper is to provide estimates leading to a direct proof 
of the exponential decay of the n-point correlation functions for certain 
unbounded models of Kac type. The methods are based on estimating 
higher order derivatives of the solution of the Witten Laplacian equation 
on one forms associated with the hamiltonian of the system. We also 
provide a formula for the Taylor coefficients of the pressure that is suitable 
for a direct proof the analyticity. 

1 Introduction 

In recent publications [66] we have given a generalization to the higher dimen- 
sional case of the exponential decay of the two-point correlation functions for 
models of Kac type. In this paper, we shall establish a weak exponential decay 
of the n-point correlation functions, and provide an exact formula suitable for 
a direct proof the analyticity of the pressure. 

Let A be a finite subset of Z d , and consider a Hamiltonian $ of the phase 
space 1R A . We shall focus on the case where $ = $a is given by 

2 

*a(x) = y + (!) 

under suitable assumptions on 

Recall that if (/) denote the mean value of / with respect to the Gibbs 
measure 

the covariance of two functions g and h is defined by 

cov(g,h) = ((g-(g))(h-(h})}. (2) 
If one wants to have an expression of the covariance in the form 

cov( 5 , h) = (Vh • w) i2(RniR „. e _ #dx) , (3) 
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for a suitable vector field w we get, after observing that V/i = V(/i — (h)), and 
integrating by parts, 

cov( 5 , h) = J(h- (/i))(V$ - V) • we-^dx. (4) 

(Here we have assumed that g and h are functions of polynomial growth ). 
This leads to the question of solving the equation 

g - (g) = (V$ - V) • w. (5) 

Now, trying to solve this above equation with w = V/, we obtain the equation 

g - (g) = (-A + V<& • V) / \ 

if) = o. J w 

The existence and smoothness of the solution of this equation were established 
in [8] (see also [66]) under certain assumptions on $. Now taking gradient on 
both sides of (6), we get 

Vg = [(-A + V$ • V) ® /d+Hess$] V/. (7) 

We then obtain the emergence of two differential operators: 

A<£ } := -A + V$ • V (8) 

and 

:= A ( £ ] ® 7d+Hess$. (9) 

Here the tensor notation means that A^ acts diagonally on the vector field 
solution to produce a system of equations. 
Thus 



cov 



(g,h) = J (A<£ rl Vg-Vh)e-*Wdx. (10) 

The operators and A^ are called the Hclffcr-Sjostrand's operators. These 
are unbounded operators acting on the weighted spaces 

L 2 (M A ,e-*rfa;) and L 2 (R A , M A , e~*G?x) 

respectively. 

The formula (10) was introduced by Helffer and Sjostrand and in some sense 
is a generalization of Brascamp-Lieb inequality as already pointed out in [1] . 
The unitary transformation 

f7$ : L 2 (M A ) — > L 2 (M. A ,e~®dx) 



will allow us to work with the unweighted spaces L 2 (R A ) and L 2 (M A ,R A ) by 
converting the operators and into equivalent operators 

mi V$| 2 A$ 
Wf = -A+^ (11) 

and 

w£ ; = I -A+J-^-l — 1 ® I + Hess*. (12) 

respectively. 

The equivalence can be seen by observing that 

w£>=e-*/ 2 o4>oe*/ 2 . (13) 

The operators wj^ and wj^ are unbounded operators acting on 

L 2 (M A ) and L 2 (R A ,R A ) 

respectively. These are in fact, the euclidean versions of the Laplacians on zero 
and one forms respectively, already introduced by E. Witten [18] in the context 
Morse theory. 

The equivalence between the operators A^ and Witten's Laplacians was 
first observed by J. Sjostrand [13] in 1996. 

2 Higher Order Exponential Estimates 

We shall consider a Hamiltonian of the form 

*(ar) = $ A (x) = y + *(&), x G R A . 

where 

\d a W\<C a , VaeN' A l. (14) 

g will denote a smooth function on R r with lattice support S g = T (^ A) . We 
shall identify g with g defined on R A and shall assume that 

\d a Vg\<C a VaGN |r| . (15) 

As in [66] , we shall momentarily assume that W is compactly supported in R A 
and g is compactly supported in R r but these assumptions will be relaxed later 
on. 

Let M be the diagonal matrix 

M = (S ijP (i)) iJeA 
where p is a weight function on A satisfying 

e -\ < < gA if • ^ • f gomc A Q QgN 

" P(j) ~ 
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Assume also that for every M as above, there exists 5 € (0, 1) such that 

(Af- 1 Hess$(a;)Ma, a) > S a a 2 , Vx E K A , Va e M A . (17) 
For instance, the d— dimensional nearest neighbor Kac model 



^a(x) = y~ 2 z2 lncosh y 2 ^ + ^ 



satisfies this assumption for v small enough. See [66] for details. 
The following theorem has been proved in [66]: 

Theorem 1 (A. Lo [66]) Let g be a smooth function with compact support on 
M. r satisfying 

\d a Vg\ <C a Va G N |r| (18) 
and is as above. If f is the unique C°° — solution of the equation 

-A/ + V$ • V/ = g - (g) 

</>L^)=0, 

then 

ieA 

n and C are positive constants. C could possibly depend on the size of the 
support of g but does not depend on A and f. 

We now propose to generalize this theorem to higher order derivatives. 

Proposition 2 If in addition to the assumptions of theorem 1, $ satisfies the 
following growth condition: for k > as above, 

*l j x il ...x th {B)e™ ({il ' m "' ik} ' Sg) < °k V.t e K A , for k > 2 (19) 

j,h,...,i k eA. 

for some Ck > not dependent on A and f, then for any k>l,f satisfies 

V fl x . ( x ) e 2K<i({n,..,* fc },S s ) < c \JxER A (20) 

' S 1 1 z k 

ii,...,ife£A 

where Ck, g > is a constant that depends on the size of the support of g but not 
on A and f. 

Proof. 

The case k = 1 being theorem 1, we assume for induction that the result is 
true when k is replaced by k < k with k > 2. 
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For k > 2 (see [8] for details), we have 

V k gM® -®tk) = (V$- V- A)^V fe /,ti®...®i fc 

fe 

+ ^ (v k f 1 t 1 ... ®Hess$ij ...<g)ife 



+ 51 <*A(ax)v$,t B (5 x )v/>. 

AUS={l,...,fe},AnS=0 
#B<fc-2 

In the right hand side of this last above equality, we have used the notation 
tj(d x )u := ( V #J u, ti ® ... ® i #J \ . 



Now fix « 2 , ...,«fe G A. Because V /(a;) — > as \x\ — > oo (see [66]), wc consider 
x a E M A that maximizes 

x < — >^2fx il ...x ik P 2 ( i i,-,ik) 

ii 

where 

p(i 1 ,...,i k ) = e Kd ^-' i ^ s °l 

Observe here that x could possibly depend on i 2 , ik £ A. 
Choose 

h = (p(h, -,ik)fx il ...x i . (x )) 

and 

t j =e ij if j = 2,...,k 
Let Mi be the diagonal matrix 

Mi = (5 sil p(ii,...,u)) sil 

and 

M,=I if j^l (21) 

in particular, we have 

V fe #,Miti ®...0M fe f fe 



(V* • V - A) ( V K /, Miti ... ® M fc ^ 



+ ^ / V fc /, Miii ® ... ® Hess^Mj-ij ® ... ® M fc i fe 



+ J! (iMA(^)V$,fMs(9x)V/) 



AuB={l fc}.An_B= 

#B<fc-2 



t M A(d x )u := (V* A f, M x t iH <g> ... ® M #A t l]#A J, Jt G A. 
As in [66], the function 

x i — ► ( V fc /(x), Miii <8> ... <8> M k t k ' 



achieves its maximum at x a . Using the notation & XiA = ^x 4< ...x 4< if A 
{£i, ...£ r } C {1, ...k} , we therefore have 

9x il ...x ik (Xo)p(il,...,ik) 2 fx tl ...x th (Xo) 
> ^""^ ,fx i - l ...x ih ( X o)fxsXi ? ...x ih (x )p(il,...,ik) ^x B x il (x ) 

fc 

+ XXX^i" ^ ...x ik (Xo)fx il ...x ik (Xo)p(il,-,ik) 2 $x s x ij (Xo) 

i=2ueA se A ^ 
+ X] ( £ v *x, A /x ll ...x llb (!Co)p(*i,-,tfc) 2 ,V/ XiB (a; ) 

AuB={l fe}.An_B=0 \ii 6 A / 

#B<fc-2 

+ ^ (v$ XiA ,^2Vf XiB (x )f Xii ... Xik (x )p(ii,...,i k ) 2 

AuB={l k}.AnB=$ \ »i£A / 

#B<fc-2 

Equivalently 

X! 9x^...x ik {x )p(ii, ■■■,ik) 2 fx tl ...x tk (x ) 

»i€A 

-XX /^ 1 ...^ fc (a;o)/a ;s ... a;ifc (a;o)p(ji,--,«fc) 2 $ 2;s a ;il (a;o) 

s£A ii£A 
fe 

+ X X X-^n-- x * ...x ife (a;o)/x il ...x ife (a;o)p(ii,-,«fc) 2 ^x s x i:J .(a; ) 

+ X X) X ®x iA xJx il ...x ik (x )p(ii,-,ik) 2 fx iB x s ( x o) 

AUB={l,...,fc},AnB=0 seA ii£A 
#B<fc-2 
l€A 

+ X X 5Z $ x iA x./x 4B x.(a; )/x 4l ...x ijl (!Co)p(*l,-,tfc) 2 

AUB={l,...,fc},AnB=0ti€A seA 
#B<fc-2 



Now taking summation over i 2 , ...,ik, we get 



X X 9x il ...xi k (x )p(ii, -,ik) 2 fx il ...x ik (x ) 

i2,...,ii;£A ii 6 A 

- S S S /^i-^ fe ( a; o)/x s -x ifc (a;o)p(ii,-,ife) 2 *x 5 x il (a; ) 

i 2 ,...,»n,€A seA ii£A 
fe 

+ X X] X] z s ...x ik (Xo)fx il ...x ik (Xo)p(il,-,ik) 2 ®x e x ij (Xo) 

i 2 ,...a k eAj=2i 1 £AseA '^Z' 

+ X X ^2^2$x iA xJx il ...x ik (Xo)p(il,-,ik) 2 fx iB x s (Xo) 

i 2 ,....ifc£A AUB={1 fc}.AnB=0 seA iiEA 

#B<fc-2 

+ X X X E $ ^ i '^B*.w/*«i-^( i °w ii '*) 2 ' 

i 2 ,...,i fc eA AUS={l,...,fe},AnS=0 *i£A seA 
#B<fc-2 
lEB 

Next, we propose to estimate each term of the right hand side of this above 
inequality. 

X XX! /x il ...x ife (a;o)/x s ...x ife (a;o)p(ii,-,ife) 2 ^x s x il (a;o) 

i 2 ,...,i fc eA sGAiiGA 

X (v/x la ...x lfc (!Co),Hes8$Miti^ 

i 2 ,...,i fe eA 

^ (M 1 V/ Xi2 ... Xife (a; ),M 1 - 1 He S s$M 1 i 1 ) 

i 2 ,...,ifc£A 

^ (*i,Mf ^ess^Miii) 

i 2 ,...,ifc£A 
i 2 ,...,i k eA 

= S ° fx il ...x ik (xo)p(iu-,ikf- 

»i,...»fc€A 



Similarly, it is easy to see that 



fe 

X EEE^i- « s ...x^^o)^!..^ (a; )p(«i,-,«fe) 2 *x 8 x i .(a;o) 

i 2 ,...,i fe eA j=2 tiEAseA V "7th"' 

> (fe-l)5 fx il ...x ik (xo)p(ii,-,ik) 2 

ii,...ifcgA 



To estimate the last two terms, we have 



E E E E \^x iA Xsfx H ...x ik {Xo)p{il,-,ikffx iB x s {Xo) 



i 2 ,...,i k EA AUB={l,...,k},AnB=0heA sGA 
#B<fc-2 



< 



E fx tl ...x tk (x )p(h,--,ik) 2 

i!,...,i k eA 



1\2 



E 

ii,...,i fe eA 



x s \%o ) 



AuB={l,...,k},AnB=9 sGA 
#B<fc-2 
ISA 



1\2 



To estimate the second factor of the right hand side of this last above inequality, 
we have 



E 

»i,...,ueA 



E E \^ iA x a {x )p{h,-,ik)fx iB x 3 ( 3 



AuB={l,...,k},AnB=® sGA 
#B<fc-2 
1GA 



< 



< 



c k E 



ti,...,U€A AuB={l 

#B<fc-2 

ieA 



ft E E 

»!,...,»* £A AuB = {l....,k}.AnB = 

#B<k-2 

ieA 



...,k},AnB=% VsGA / 
(E*x iA s>o)A*l. ■■■.**)) X 

VseA / 



< 



(5Z^(*l.-»*fc)/x iB x.(*o) 
VsGA / 

ft E E (E'": , ^ ') 

»i,...,i)teA AuB = {l,...,fc},AnS=0 VsGA / 



#B<k-2 

ieA 



(^ e 2 K d({i r .jeB}u{s},S g )j2 ^ ^ Xo A 
VsGA / 



< Cfc. 

This last inequality above follows from the induction assumption and that of 
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Thus, 



E E EE* 

fx il ■■■Xi k (%o)p(i>li ik) fxi B x s { x o) 



,**€A AUS={l,...,fc},AnS=0iieAseA 
#B<fc-2 



- - E fl il ...x ik {x )p{ii,-,ikf 

ii,...,»i,eA 



1 1\2 



Similarly, we have 



E E E E^a^^b^^ )^!-^^ )^ 1 '-'^) 2 



i2,.")*fc£A AUB={l,...,fe},AnB=0 ii£A seA 
#B<fc-2 
1GB 



> -c fc 

We then finally get 



E fx zl ...x Zk ( x o)p(h,--,ikf 

ii,...,ifc€A 



1\2 



X] 9x H ...x ik (xo)p(h, -,ik) 2 fx il ...x ik (x ) 

ii,...,i k eA 

> k5„ fl il ...x ik { x o)p{h,...,ikf 



ii,...i k £A 



E fx tl ...X> k (Xo)p(h,-;ik) 2 

ii,...,i k £A 



1\2 



If 



E /x 4l ...x ij ,(*o)P(*l.-.*fc) 2 = 
ii ,...,ifc GA 

then there is nothing to prove, otherwise we have, after using Cauchy- Schwartz 
and dividing by 

E fx tl ...x tk (xo)p(h,-,ik) 2 , 

ii,...,i k eA 



< 



1/2 



E fx tl -x tk ( x o)p(h,-,ik) 7 

\ii,...ik&A 

^ E 9l H ... Xik {xo)\ +C Kg 
\h,...,i k eA I 



< Ck,g- 
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3 Relaxing the Assumptions of Compact sup- 
port 

As in [8] , we consider the family cutoff functions 

X = Xe (22) 
(e G [0, 1]) in C~(R) with value in [0, 1] such that 

X = 1 for |t| < e- 1 

\x {k) (t)\ < C k -^ for fee N. 

We then introduce 

* e (a;) = x e (kl)* (*) (23) 

and 

Je (i)= Xe (|i|) fl (i) xeR r . (24) 
A straightforward computation (see [66]) shows that ^f s (x) and g E (x). satisfy 

|9 Q V* e | <C a + Q , A ( £ ), Va G N' A ' . (25) 

and 

\d a Vg e \ <C a + O aA (s), \ta G N^, (26) 

and that 

M- l H.ess$ e (x)M > 5', < 5' < 1. in the sense of (17) (27) 

It then only remains to check that 

V tf ? ( a .) e 2i«i({ii,...,i»},s») < Ck + O k A (e) G K A , Vfc > 2 

j,*i,...,»fceA 

(28) 

where Cfc is a positive constant that does not depend on / and A. 

* e (ar) = X e (r)tf(a;) 
Let a be such that |a| > 3. Using Leibniz's formula, we have 



(29) 



< |0 Q x £ (r)*| + + E |^x £ (r)a Q -^| . (30) 



/3<a 

' a 



Assuming that \I/(0) =0 and write 



<&(x) = f x ■ V^(sx)ds 
Jo 
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Now using the fact that 
we have 



ji £A- 

< C\rd<*Xe{r)\. 
\d a x £ (r)V(x)\ =O aA (e). 



Finally, using the fact that 

d^xAr) = O {e) for every > 1, 
it is then easy to see that 



(31) 



O aA {e). 



(32) 



Thus 



V tf* ,y K d({h,..., i k },s g ) < Ckq + O kA (e) Va;eR A ,Vfc>2. 

],n,...,i k eA 

(33) 

Now using the arguments developed in [8] (see also [66]) about the conver- 
gence of the corresponding solutions as e — ► 0, we obtain: 

Proposition 3 If g(0) = X P(0) = 0, then Proposition 2 holds without the as- 
sumptions of compact support on ^ and g. 

4 The n-Point Correlation Functions 

The higher order correlation is defined as 

(9u -,9k) ■= - (9i))-(9k - (9k)))- (34) 

For simplicity we shall take k = 3 and $ is as in proposition 2. 

Let gi,g2, and g$ be smooth functions satisfying (15) and fi i = 1,2,3 shall 
denote the unique solution of the system 



-Afi + V$ ■ V/i = 9i - (g % ) 



(fi) 



0. 



1.2 („) 



(35) 



Recall that 



v/ i = ^ 1) " 1 v ft . 
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For an arbitrary smooth function c, it is easy to see that 

{c(x) (ft - <ft») = <V/i ■ Vc) . 

A direct computation shows that 

(9i, 92,9s) = (V/ 3 • (Hess/i) V.g 2 ) + (V/ 3 • (Hess.g 2 ) V/i> 

+ (V/ 2 • (Hess/i) V.g 3 > + (V/ 2 • (Hess 53 ) V/i) . 

Let us now estimate each term of the right and side of this equality. 
Using Cauchy-Schwartz, and proposition 2, it is easy to see that 

|(V/ 3 • (Hess/0 V<fc>| < Ce- Kld ( 5 ^i) 

|(V/ 3 • (Hess 52 ) VA) | < Ce- Kld (^' s "), 
|(V/ 2 • (Hess/0 V 53 )| < C e - Kl < s »3' s '0 

and 

|(V/ 2 • (Hess 53 ) V/i)| < Ce- Kl < s ^ s ») 

Here the constants C only depends on the size of the support of the gi's. and 
Ki > 0. 
Thus 



\(9i,92,93) \ < c 



-K 1 d(Sg 2 ,Sg 1 ) _|_ e -Kld(Sg 3 ,S gi ) 



If gi = Xi, gi = Xj, and g 3 = Xk, we obtain 

\((xi - (xi)) {xj - ( Xj )) (x k - (x k )))\ < C 



e -md(i,j) _|_ e -md(i,k) 



Thus if d > 1, we obtain this weak exponential decay of the truncated correla- 
tions in the sense that the exponential decay occurs as you simultaneously pull 
the spins away from a fixed one. Note that in the one dimensional case, we 
obtain a stronger exponential decay due to the fact that 

i < j <k => d(i, k) = d(i,j) + d(j, k). 

This was already pointed out in [8]. 



5 The Analyticity of the Pressure 

In this section, we attempt to study a direct method for the analyticity of the 
pressure for certain classical convex unbounded spin systems. It is central in 
Statistical Mechanics to study the differentiability or even the analyticity of the 
pressure with respect to some distinguished thermodynamic parameters such as 
temperature, chemical potential or external field. In fact the analytic behavior of 
the pressure is the classical thermodynamic indicator for the absence or existence 
of phase transition. The most famous result on the analyticity of the pressure 
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is the circle theorem of Lee and Yang [28] . This theorem asserts the following: 
consider a {—1, 1} —valued spin system with ferromagnetic pair interaction and 
external field h and regard the quantity z — e h as a complex parameter, then all 
zeroes of all partition functions (with free boundary condition), considered as 
functions of z lie in the complex unit circle. This theorem readily implies that 
the pressure is an analytic function of h in the region h > and h < 0. Heilmann 
[29] showed that the assumption of pair interaction is necessary. A transparent 
approach to the circle theorem was found by Asano [30] and developed further by 
Ruellc [31], [32], Slawny [33], and Gruber et al [34]. Griffiths [35] and Griffiths- 
Simon [36] found a method of extending the Lee- Yang theorem to real-valued 
spin systems with a particular type of a priory measure. Newman [37] proved 
the Lee- Yang theorem for every a priory measure which satisfies this theorem in 
the particular case of no interaction. Dunlop [38], [39] studied the zeroes of the 
partition functions for the plane rotor model. A general Lee- Yang theorem for 
multicomponent systems was finally proved by Lieb and Sokal [40]. For further 
references see Glimm and Jaffe [41]. 

The Lee- Yang theorem and its variants depend on the ferromagnetic char- 
acter of the interaction. There are various other way of proving the infinite 
differentiability or the analyticity of the pressure for (ferromagnetic and non 
ferromagnetic) systems at high temperatures, or at low temperatures, or at 
large external fields. Most of these take advantage of a sufficiently rapid decay 
of correlations and /or cluster expansion methods. Here is a small sample of rele- 
vant references. Bricmont, Lebowitz and Pfister [42], Dobroshin [43], Dobroshin 
and Sholsman [44], [45], Duncau et al [46] , [47] , [48] , Glimm and Jaffe [41], [49], Is- 
rael [50], Kotecky and Prciss [51], Kunz [52], Lebowitz [53], [54], Malyshev [55], 
Malychev and Milnos [56] and Prakash [57]. M. Kac and J.M. Luttinger [58] 
obtained a formula for the pressure in terms of irreducible distribution functions. 

We propose a new way of analyzing the analyticity of the pressure for certain 
unbounded models through a representation by means of the Witten Laplacians 
of the coefficients in the Taylor series expansion. The methods known up to 
now rely on complicated indirect arguments. 

6 Towards the analyticity of the Pressure 

Let A be a finite domain in TL d (d > 1) and consider the Hamiltonian of the 
phase space given by, 

&(x) = $a(x) = — + tf (a;), x e K A . (36) 

where 

\d a Vt>\<C a , VaeN |A| , (37) 

Hess$(a;) > <5 Q , < S Q < 1. (38) 

Let g is a smooth function on M r with lattice support S g — T. We identified 
with g defined on R A by 

g(x) = g(x r ) where x = (xi) ieA and x r = {x l ) ier (39) 
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and satisfying 

\d a Vg\ <C a Vae N |r| (40) 

Under the additional assumptions that is compactly supported in R A and g is 
compactly supported in R r , it was proved in [66] (see also [8]) that the equation 

f -A/ + VS-V/ = <?-<<?) 
I </W) = 

has a unique smooth solution satisfying V fc /(x) — > as |x| — > oo for every 
k > 1. 

Recall also that V/ is a solution of the system 

(- A + V$ • V) V/ + Hess$V.f = Vg in R A . (41) 

As in [66] and [8], these assumptions will be relaxed later on. 
Let 

*i(x) = *(x)-tg(x), (42) 
where x — (xi) i£ A, and assume additionally that g satisfies 

Hessg < C. (43) 

We consider the following perturbation 

6 A (t) = log J dxe~**W 

Denote by 



(44) 



and 



Z t = J dxe-*^ (45) 
[■dxe-*^ 

< ■ >t,A= 1 y • (46) 

7 Parameter Dependency of the Solution 

From the assumptions made on $ and g, it is easy to see that there exists T > 
such that or every t £ [0, T), ^(x) satisfies all the assumptions required for the 
solvability, regularity and asymptotic behavior of the solution f(t) associated 
with the potential &\(x). Thus, each t € [0, T) is associated with a unique 
C°°— solution, f(t) of the equation 

r ^/(*)=5-( ff u, 

1 (/(*)>^ w = o. 



Hence, 

A$v{t)=Vg (47) 
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where v(i) = V/(t). Notice that the map 

t .— > v(t) 

is well defined and 

{v(t):t€[0,T)} 
is a family of smooth solutions on R A satisfying 

d a v(t) -> as \x\ -> oo Va£ N |A| and for each t e [0, T) 

and corresponding to the family of potential 

{*i:te[0,T)}. (48) 

Let us now verify that v is a smooth function of i G (0, T).We need to prove 
that for each t € (0,T), the limit 

hm V(t + £) - vW 

e^O £ 

exists. Let 

v£(t) = v(* + e)-v(*), 

e 

We use a technique based on regularity estimates to get a uniform control of 
v e (t) with respect to e. 

With e small enough, we have 



= -A 



v(t + e) - v(t) 



St, 



^_ V$*+ £ ■ Vv(t + g) - y$« ■ Vv(f) 



+ 

Equivalently, 

-A 



Hess$ t+e v(t + e) - Hess$'v(i) 



v(t + e) - v(t) 



+Hess$ t+e 



£ 

v(t + £) - v(t) 

£ 



V$ t+£ ■ V [v(t + g) - w{t)} 



and 



= -( jv(*)-( )-Vv«/, 

= Hess.gv(i) + Vg • Vv(t) 



-Av e (t) + V$* +£ • Vv e (i) + Hess$* +£ v e (t) 
Hess.gv(t) + • Vv(i) 
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Let w(t) be the unique C°°— solution of the system 

- Aw(t) + V$* • Vw(t) + Hess$*w(t) = Hess.gv(t) + Vg • Vv(t). (49) 
Recall that the unitary transformation C/$t+ Ej allows us to reduce 



-Av £ (t) + V$* +e • Vv £ (f) + Hess$*+ £ v £ (i) 
= Hess.9v(t) + V# • Vv(i) 



(50) 



into 

/ lV$ t+,r l 2 A$ i+£ \ 

- A+ J 1 V £ + Hess$ t+e V £ = , N 

V 4 2 I (51) 



[Hessflv(t) + V.g • Vv(i)] e"* / 2 
where V £ = v £ (t) e -* t+e / 2 . 

Remark 4 TTiis unitary transformation already mentioned in the introduction 
was introduced in the proof of the existence of solution (see [66] ) to avoid work- 
ing with the weighted spaces L 2 (R A , R A , e~^dx). The proof was based on Hilbert 
space method. The method consists of determining an appropriate function space 
and an operator which is a natural realization of the problem. In this particular 
problem, the function spaces to be considered are the Sobolev spaces _B|(M A ) 
defined by 

B|(R A ) = {ue L 2 (R A ) : Z^d a u E i 2 (M A ) V I + |a| < k) . 

where 

|V$ 

(52) 



2 

These are subspaces of the well known Sobolev spaces W k ' 2 (R A ) 7 k £ N. 



Taking scalar product with V £ on both sides of (51), we get 

+ / Hess$*+ £ V £ • V e dx = 



(53) 

/ [Hessgv(t) + Vg • Vv(t)] e -** +E /2 . V £ dx. 

Now using the uniform strict convexity on the left hand side and Cauchy- 
Schwartz on the right hand side, we obtain 

||V £ || B o< C t for small enough e. (54) 

We then deduce that 

-A+ 1 — A V £ = q £ (55) 
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where 

A<J>*+ £ 

q £ = [Hess.gv(t) + V.g • Vv(t)] e~* /2 + V £ - Hess$ t+£ V £ (56) 

is bounded in _B° uniformly with respect to e for e small enough. 

Taking again scalar product with V £ on both sides of (55) and integrating 
by parts, we obtain 



|vv £ || 2 2 



|V$ t+£ | 



<ll?s|| i2 ||Vl i2 (57) 

L 2 



It follows that V e is uniformly bounded with respect to e in B\ t+E for e small 
enough. 

Next, observe that 

( A+ \J^L) V £ = ~ £ (5g) 



where 



Qe = ? £ - J ^ ^V £ + ^ >- V £ (59) 

= , £ -£l^v £ -£^^v £ (60) 

is uniformly bounded in B° with respect to e for small enough e. Using regu- 
larity, it follows that for small enough e, V £ is uniformly bounded in B^ t with 
respect to e.This implies that q e is uniformly bounded in B\ t for e small enough. 
Again, we can continue by a bootstrap argument to consequently get that for e 
small enough, V £ is uniformly bounded in B^ t with respect to e for any k. 
Let 

V = w(t)e-*V2. 

We have 

/ |V$ 4 | 2 A$*\ 

Jv + H«b*«V (61) 

= [Hessgv(t) + Vg • Vv(i)] e"*'/ 2 
Now combining this equation with (51), we obtain 



-A+^-^ — I (V £ - V) + Hess** (V £ - V) 



= - [Hess.gv(i) + V.g • Vv(t)] e"**/ 2 + [Hess<?v(i) + Vg • Vv(t)] e -* t+e / 2 
+ 



|V$*r |V$* +£ | 2 \ /A$* A$ t+£ 



4 / V 2 2 

+ (Hess$* - Hess$ t+£ ) V £ . 

(62) 
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Now let us check that for small enough e, the right hand side of (62) is 0(e) 
in B°. 

For the first term, we have 

- [Hessgv(t) + V.g • Vv(t)] e"**/ 2 + [Hessgv(t) + V.g • Vv(t)] e -* t+E /2 
= [Hess.gv(t) + V.g • Vv(t)] e"**/ 2 ( e e 9 /2 _ ^ 
~ § [Hess 5 v(f) + V.g • Vv(t)] ffe"**/ 2 

Thus for e small enough 

[Hessgv(t) + V ff • Vv(t)] e"**/ 2 + [Hess ff v(t) + Vg • Vv(t)] e -* t+E / 2 
For the second term, we have 

|V$ 4 | 2 |V$ t+£ | 2 



< 



4 4 
= I (|v$*| + |V$* +£ |) (|V$*| + |V$* +£ |) 

< | \Vg\ (|V$*| + |V$* +£ |) 

< ||V ff |(2|V$*|+ £ |V.g|) 

Using now the fact that V e is uniformly bounded in B^ t with respect to e 
for any k, we see that the second term of the right hand side of (62) is 0(e) 
in B^ t . The last two terms of the right hand side of (62) are obviously 0(e) in 

From the same regularity argument as above, we get that V £ — V is 0(e) in 
i?| t . Again iterating the regularity argument, we obtain that for small enough 
e, V £ — V is 0(e) in B\ t for every k. We have proved: 

Proposition 5 Under the above assumptions on <I> and g, there exists T > 
so that for each t 6 (0, T), v e (t) converges to w(t) in C°° . 

Remark 6 The proposition establishes thatv(t) is differentiable int and 
is given by the unique C°° — solution w(t) of the system 

- Aw(t) + V$* • Vw(t) + Hess$*w(t) = Hess.gv(t) - Vg • Vv(t). (63) 

Iterating this argument, we easily get that, v(t) is smooth in t G (0,T). 

Now we are ready for the following: 

8 A Formula for the Taylor Coefficients 

First observe that for an arbitrary suitable function f{t) = f(t, w) 

g- t < f(t) >t,A=< / '(*) >t,A +cov(/, g). (64) 
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Hence, 

Let 
Thus, 



-g t < f{t) >t, A =< / '(*) >t, A + < 4V" (V/) • V.g > t , A . 



^/ == 4V" 1 (V/) • V.g. 



at 

d 



< f{t) >t,A=< 



at 



A? / >t,A • 



(65) 
(66) 
(67) 



The linear operator — + A g will be denoted by H g . 

To obtain a formula for the coefficients in the Taylor expansion of 



/ 



dxe-**™ 



e A (t) = log 

we first the derivatives of 6\(t) in terms of H g 

( d \° 

9' A (t) =< g > t ,A=< (g- t + A 9j 9 >t,A=< H g9 >t,A] 

Kit) = §1<9 >t,A=< 4 1 *' 1 ( V -9) ' V -9 >t,A=< + Ag] g > tlA ; 



(68) 



St 



<?a(*) = ^<4 i r 1 (v fl )-v ff > t ,A=<^(4 i r 1 (v ff ).v 5 )>t,/ 
+ < (4 i r 1 v(4r 1 (v fl )-v 5 ))-v fl > tlA 



By induction it is easy to see that 
d x 



*£°(t) =< ( m + A g ) g >t, A =< >t,A (Vn > 1) 



Next, we propose to find a simpler formula for 9^\t) that only involves A g . 

H g g = 4 1 ?" 1 (V 5 ) • V 5 



H 2 g g - |v/ • V.g + (a^ V (4V"' (V.g) • V.g)) • Vg 



(69) 



where / satisfies the equation 



V/ = 4V" 1 iVg) . 



(70) 
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With v(t) = V/, as before, we get 



^ V/ • Vg = A%r (Hessgv(t) + Vg ■ Vv(i)) • Vg 



and Hg becomes 

U gS - A <5>t 



(Hess 3 v(t) + Vg ■ Vv(i)) + V (i™"' (V$) • V 5 )] • V.9 



AW _1 2V(^)-V ff 



= 2^. 9 . 
Proposition 7 7/ 



Mt) = iog 



-where 



$*(x) = $a(x) - tg(x) 



is as above then 0^\t), the nth— derivative of 9 \(t) is given by the formula 



and for n > 1 



6' A (t) =< g > t , A , 



9^(t) = (n-l)\<A n g- 1 g> 



t,A 



Proof. We have already established that 



e^W^H^g > tA forn>\. 

It then only remains to prove that 

H n g - x g = (n-\)\A n g- x g forn>l. 

The result is already established above for n = 1,2,3,. By induction, assume 
that 



m- 1 g=(n-l)\A n - 1 g. 



if n is replaced by n < n. 



H v9 = [^+A g )((n-iy.A-- 1 g) 



Now 



A V9 = 



A^rV(A^g) 



V.9 
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where 



We obtain, 



A 



4>* 



n-2. 



at 

We then have 

— A n - X n - 
dt 3 9 



£ V^„ = (|v^- 2 .g + Hess.gV^„ + V.g ■ V (V^; 



Thus, 



Vy„ • V.g 



, ; g + Hess.gV^„ + V.g • V (V<^„) 

V.g 



V.9 



4 1 ?" 1 (|v^- 2 .9 + V(V^-V 5 ) 



= A„ 



r 2 



1 



.4 H"("- 1 ) 



(n-2)! 9 9 
^-l^A g ((n-l)U r ^) 
(n-l)A^. 



(from the induction hypothesis) 



(still by the induction hypothesis) 



H n g g = (n-l)!(n-l + l)AJ ff 
= n\A n g g 



Proposition 8 If g(0) = 0, t/ien the formula 

e^\t) = {n-l)\<A n g - 1 g> t ^ n>2 
still holds if we no longer require ^ and g to be compactly supported in M A . 
Proof. As in [8] , consider the family cutoff functions 

X = Xe (71) 

(e £ [0, 1]) in C^°(M) with value in [0, 1] such that 

X = 1 for |i| < e- 1 



for k E 
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We could take for instance 

Xs (t) = m*\t\) 

for a suitable /. 

We then introduce 

* c (a;) - Xe(kl)*, x E R A (72) 

and 

g e (x) = Xe (\x\)g x e R r (73) 

One can check that both \I/ £ (x) and g £ (x) satisfies the assumptions made above 
on \t and g. Now consider the equation 

- Af £ + V$* • V/ e = g e - < g e > t , A . (74) 

which implies 

(- A + V$* • V) ® v £ + Hess$* v £ = V.g e (75) 

where 

It was proved in [8] that v e = A^ t ' V.g e converges in C°° to A^ t ' Vg as 

£ -> 0. ■ 

Proposition 9 Let 

fee i/ie finite volume Pressure. 

Denote by a n (n > 2) the nth Taylor coefficient. We have 



a n = 



< A n - X g > A 
n|A| 



Remark 10 This formula for a n gives a direction towards proving the analytic- 
ity of the pressure in the thermodynamic limit. In fact one only needs to provide 
a suitable C n estimate for < A™~ x g >\ . 

9 Some Consequences of the Formula for nth— Derivative 
of the Pressure. 

In the following, we shall additionally assume that 

Vg(0) = 0, and 
V$^(0) = for alii e [0,T). 
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When n = 1, we recall that A° g g = g, 

6' A (t) =< g > t . A 

and if 

v(t) = V/ = A$~ Vg, 

then we have 

(-A + • V) ® v(t) + Hess$^v(t) = Vg 

Again the tensor notation means that (—A + V$ A • V) acts diagonally on the 
components of v(i). 

As in [8] v(t) is a solution of the equation 

9 =< 9 >t.A +v(t) • - divv{t). (76) 

Using the assumptions above, we have 

#a(*) = < 9 >t,A 
= diuv(i)(0). 



Similarly, the formula 

implies that 
where 



»W(t) = (n - 1)! < >,,, 

^ n) (t) = (n-l)!dwv n (t)(0), 
v n (i)=4 1 t ) " 1 V(A"- 1 5 ). 
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